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Abstract 

We present new intuition behind Grover's quantum search algorithm by means of a Hamil- 
tonian. Given a black-box Boolean function / : {0, 1}" {0, 1} such that f{w) = 1 for exactly 

■ one w £ {0, ll", Grover [Q describes a quantum algorithm that finds w in 0(2"/2) time. Farhi 
I & Gutmann |^ show that w can also be found in the same amount time by letting the quantum 

■ system evolve according to a simple Hamiltonian depending only on /. Their system evolves 
I along a path far from that taken by Grover's original algorithm, however. The current paper 

presents an equally simple Hamiltonian matching Grover's algorithm step for step. The new 
Hamiltonian is similar in appearance from that of Farhi & Gutmann, but has some important dif- 
ferences, and provides new intuition for Grover's algorithm itself. This intuition both contrasts 
with and supplements other explanations of Grover's algorithm as a rotation in two dimensions, 
and suggests that the Hamiltonian-based approach to quantum algorithms can provide a useful 
' heuristic for discovering new quantum algorithms. 



1 Introduction 

Quantum algorithms can, in theory at least, solve useful problems faster than classical algorithms. 
Two primary families of quantum algorithms in this regard are algorithms for factoring and discrete 
log and Grover's search algorithms with quadratic speed-up |Q, |5|. 

There are many variations on Grover's original algorithm — counting, starting with partial data, 
multiple targets, et cetera. The algorithm is also surprisingly robust; although the original algo- 
rithm uses the Walsh-Hadamard transform, essentially any unitary operator will do just as well []. 
Starting with a simple condition on what transform is used, we will show how Grover's algorithm 
arises from a particularly simple — almost naive — intuition about quantum algorithms. Our ideas 
also generalize to variants of Grover's algorithm. 

There have been good explanations in the literature [^, ^ of how and why fast quantum search 
works: the initial state is slowly rotated (in two complex dimensions) into the target state by 
repeatedly applying a special operator known as Grover's iterate. Using a more physics-based 
approach, Farhi Sz Gutmann |§ describe an "analog" version of quantum search by means of a 
simple, time- independent Hamiltonian which transforms any initial state \a) into some prespecified 
target state \w) in optimal time, provided that |o") and \w) are not orthogonal. Their analog 
algorithm rotates the initial state into the target state in the same time that Grover's "digital" 

* Supported in part by NSF Grant CCR-9996310. Email feimer@cs.sc.edu. 



1 



algorithm does, yet their rotation strays far from the intermediate states reached in the original 
algorithm by applying Grover's iterate. 

We show here that a simple, time-independent Hamiltonian for a system of qubits results in 
time evolution matching Grover's iterate exactly. This Hamiltonian also provides a nice, simple 
insight into the workings of the algorithm that is new, to the best of our knowledge. 

Our Hamiltonian bears some resemblance to that of Farhi &: Gutmann, although ours was 
conceived independently. Ours differs from theirs in important respects, however, and may not be 
as plausible physically, but it does closely coincide with the iterations of Grover's algorithm, and 
thus gives a much closer simulation of a digital quantum circuit by an analog process and vice 
versa. While Farhi & Gutmann's Hamiltonian is appealing from a physical point of view, ours is 
appealing from an algorithmic perspective. The heuristic it gives for Grover's algorithm suggests 
that other digital quantum algorithms might be found by first looking at analog versions. (More 
work on analog algorithms has been done recently, see for example.) 

1.1 Structure of the Paper 

We give mathematical preliminaries in Section |2|, including a brief description of Grover's algorithm 
as described in In Section |3| we show how Grover's algorithm arises from a simple-minded ap- 
proach to quantum search. The operator we describe there corresponds directly to our Hamiltonian, 
which in Section |^ we compare with that of Farhi and Gutmann and show how it generates 
Grover's iterate. Most of our work was done independently of [^] before it came to our attention, so 
our approach to the problem is different. In Section |5|, we suggest a Hamiltonian-based approach 
to quantum algorithms in general, and we present open problems. 



2 Preliminaries 

We work with linear operators over a fixed A^-dimensional Hilbert space. A standard norm on 
operators is defined as 

1^1 = sup \Av\, 

\v\=l 

where A is an operator and | • | on the right hand side is the standard hermitian norm on the Hilbert 
space. This norm on operators satisfies \AB\ < \A\ ■ \B\. Clearly, all unitary operators have unit 
norm. The exponential map on operators is defined as 

e^ = / + ^ + |r + |r + ---' W 

where I is the identity operator. 

If A is skew hermitian {A^^ = —A), then is unitary. Conversely, for any unitary U there is a 
skew hermitian A such that U = e^. As in the case with the exponential function on scalars, we 
also have ^ 

e^= lim (i+t) ■ (2) 

For an n-qubit system we assume the standard basis of states \i) indexed by classical bit con- 
figurations i E {0, 1}"'. We use lower-case Roman letters to label basis states, and lower-case Greek 
letters to label other (arbitrary) states in the Hilbert space. 
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2.1 Grover's Search Algorithm 

Here we briefly review Grover's search algorithm. Fix an integer n and let N = 2^. Let / be 
a Boolean- valued function on n-bit strings such that f{w) = 1 for exactly one w (the target), 
< w < N (identifying strings with integers). A simple version of Grover's algorithm is to find w 
via a quantum algorithm where inputs to / are stored in n qubits, and / is available as a black 
box function (oracle) that can be queried by the algorithm. Alternately, we may assume that /(y) 
is efficiently computable given y, and embed the computation of / into the quantum circuit. 

In this setting, Grover's algorithm (as described in [^, |6| or Q) uses three n-qubit unitary 
transforms: 

1. an arbitrary, easy-to-compute U such that (t(;|C/|0) ^ 0, 

2. the selected inverter Iq = Eo<i<Af (-1)*"° 10(^1 = I - |0)(0|, and 

3. the selected inverter I^u = I — 2 {w) {w\ . 

(Here, the formula i = in the exponent stands for its numerical truth value — 1 for true, for 
false.) These combine to form Grover's iterate 

G = -UIqU-^I^. (3) 

By adjusting U by an appropriate phase factor, we can assume that {w\U\0) = x for some real 
X > 0. This adjustment leaves G unchanged. 

Suppose / is as above with w unique such that /(w) = 1. The algorithm starts in the state |0) 
(all qubits cleared), then U is applied to get the state 

\^P) = U\0). (4) 

Next, G is applied repeatedly to 1-0), approximately [^] times. At this point, the system will be 
very close to the state \w), so when we now measure the qubits we get w with high probability. 
Note that 

i 

SO can be simulated easily given access to / alone and some extra work qubits. 

In Grover's original presentation, U = is the Walsh-Hadamard transform on n qubits, and 

so 

i 

whence, x = {w\ip) = 2""/^, which yields the quadratic speed-up in the search. 

3 Quantum Search Revisited 

The point of this section is to show how one might stumble upon Grover's algorithm by taking a 
simplistic, almost naive, approach to quantum search. The intuition here is not geometric, as it is 
with Jozsa rather, it is purely algorithmic in fiavor. 

We start with the basic observation that if j4 is a skew hermitian operator {A^ = —A) and 
< e << 1, then I + eA approximates e*^"^, which is unitary. Therefore, I + eA approximates a 
plausible step in a quantum computation. By (^), we can approximate the action of on a state 
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by repeatedly applying / + to the state roughly 1/e times. The smaller e is, the better the 
approximation. (In general, it is not certain that e*^^ is renderable by a small quantum circuit; it 
will be in the present case, though.) 

A simple example is when A = \i){j\ — \ j){i\ for some i,j £ {0, 1}", i ^ j. Applying I + eA to 
a state \ip) = J2i 10 gives 

(/ + eA) Iv?) = |(^) + eaj \i) - eai \j) . 

The operator alters \ip) (viewed as a column vector) by adding an e fraction of its jth component 
into its ith component, and in exchange, subtracting an e fraction of its ith component from its 
jth component. In a sense, we are moving probability amplitude from state |j) to state With 
arbitrary A, this swap may take place between many pairs of components of \(p) at once. 

Suppose we are given n, A^, /, and w as in Section |1[ We start in the state = N~'^/'^ Y.i K), 
which we would like to transform to the target state \w). A promising way to do this, given our 
considerations above, is to pile positive probability amplitude onto \w) while taking it away from 
all the other states evenly. The real skew symmetric operator that does this is 
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expressed in the {|i)} basis, where the nonzero entries are all in the wih. row and wih column. We 
see that 

1 ... -e ... 



I + eA 



1 -e 
e 1 e 
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The e's on row w have the effect of giving probability amplitude to \w) while removing it from all 
the other states evenly (the column of —e's). The probability amplitude of \w) gains at the expense 
of an e fraction of all the other probability amplitudes. From this it is clear that if we start in state 
where all the probability amplitudes are equal, and apply I + eA (for some small e) the right 
number of times, eventually the state \w) will dominate. 
We note that, using bracket notation, 

A -- 



N{\w){i;\-\i;){w\). 



The operator ieA acts as a Hamiltonian for the time evolution of the system from {ip) to \w). As 



we'll see in the next section, for the right value of e, 
iterate. 



e*^ is exactly two applications of Grover's 
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4 Hamiltonians 



In this section, we give a Hamiltonian for Grover's algorithm, that is, an operator H such that 
^-tHt follows the course of the algorithm as t increases. It is clear both by geometric considerations 
and by the last section that such an operator must exist. H is analogous to a previous Hamiltonian 
H' for quantum search found by Farhi & Gutmann Q which does not match Grover's algorithm. 
We first briefly describe their results, then describe our Hamiltonian using their framework. 

4.1 Farhi & Gutmann's Hamiltonian 

We are given n, A^, / and w as above. Farhi & Gutmann describe a physical, analog way to do 
quantum search by first assuming that a Hamiltonian 

= E \ w){w\ 

is available that distinguishes the target state \w) from all others by giving it some positive energy 
E (the other basis states have energy 0). Let \a) be some arbitrary unit vector in the Hilbert space 
(the "start" state). We assume \a) is easy to prepare, so for example, \a) may be \tp) of equation 
@. The goal is to evolve from \a) into \w). To search for the state \w), we are allowed to add 
some "driver" Hamiltonian Hd to Hyj, provided that Hd does not depend on the actual value of 
w at all. They choose Ho = E\a){a\, so their Hamiltonian is 

H' = Hd + H^ = E {\a) {a\ + \w){w\) , 

where E is some arbitrary positive value in units of energy. If \a) and \w) are not orthogonal, then 
we can assume as before that {(y\w) = {w\(t) = x for some x > by adjusting \a) by an appropriate 
phase factor. 

Applying e~*^ * to the start state \a) gives the time-evolution of the system,]^ which stays in 
the two-dimensional subspace spanned by \a) and \w). Restricting our attention to this subspace, 
it is easy to see that H' has eigenvalues E[l it x) with corresponding eigenvectors 

1+') = (2 + 2x)-V2(|^) + |^)), 
I-') = {2-2x)-'lWa)-\w)). 

A straightforward calculation yields 

e"^^'* \a) = e"'^* [cos{xEt) \a) - i sm{xEt) \w)] . (5) 

When t = ■7t/{2Ex), we have 

as desired. 

Farhi &: Gutmann observe that if the unit vector \a) is chosen at random, then the expected 
value of X is N~^/'^, making t = 0{N^/'^ / E). For constant E, this time is the same order of 
magnitude as Grover's algorithm. They show that their time evolution is optimal up to an order 
of magnitude for any tD-independent driver Hamiltonian Ho, even one that varies with time. 

^The evolution of a quantum system under a time-independent Hamiltonian H is actually e"'^*''''. We choose 
units so that 7i = 1, and so Et is a unitless quantity. 
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4.2 Another Hamiltonian 



The time evolution of the system according to H' strays far from the intermediate steps Grover's 
algorithm. There surely is a Hamiltonian, however, whose time evolution matches the steps of 
Grover's algorithm exactly, since each step of Grover's algorithm essentially amounts to a rotation 
in a two-dimensional space. We show that this Hamiltonian can be described very simply: the 
operator ieA mentioned at the end of Section ^ is exactly the Hamiltonian in question, for an 
appropriate value of e which we will calculate. 

The fact that Grover's iterate can be rendered by a small quantum circuit then tells us that 
our intuition of Section ^ is justified: the incremental application oi I + eA indeed corresponds to 
a legitimate quantum algorithm. 

Given n, A^, /, \w), |cr), Hy^ and H£) as above, with {w\a) = a; > 0, we define the Hamiltonian 

H='^ [H^,Hd] = 2iEx{\w){a\ - \a){w\). 
The rest of this section is devoted to proving the following 

Theorem 1 Assume the special case where \a) = {ip) and E = 1. Restricted to the , \ w)) -plane, 
e^*^ approximates Grover's iterate G to within 0{N^^'^) in norm. In fact, e~^^^° exactly matches 
G where 

vr — 2 arccos x 
= 2x^1 -x^ ■ 

On the whole Hilbert space, e~^*^ approximates G^ to within 0(A^^^/^), and e~^*-'^*o = G^ . 

For the moment, we allow E to be any positive value and \a) an arbitrary unit vector with 
< {w\a) = x < 1. Restricting our attention to the subspace spanned by \a) and {w), and letting 
= arccos X, the eigenvalues of H are seen to be ±^E sin29 with corresponding eigenvectors 



l±) = ( e^'' \a) - \w] 



1 

A/2sm( 

Setting 7] = E sin 29 = 2Exsm9, a routine calculation shows that 

e-*^*|cT) = [sm{9 - r]t) \a) + sm{r]t) \w)] , (7) 

sin^ 

e~'"^\w) = [- sm{r]t) \a) + sm{9 + r]t) \w)] . (8) 

sin^ 

If X is smah, 9 will be close to 7r/2. For t = 9/r] = 9/{Esm29) = tt/{2Ex) we have e~*^* \a) = \w). 
That is, the system finds the target state in roughly the same time as with H' . 

Comparing (||) and (Q), we see that the quantum system evolves significantly differently under 
the two Hamiltonians H' and H — by more than just a global phase factor. We now show how the 
latter evolution, run for a short time interval, matches a single step of Grover's algorithm (one 
application of G). We now assume E = 1 and \a) = \Tp) = U \0) given by equation (Q), with G 
given by (^). We again set x = {cr\w) = {^jIw) = cos 6* > 0, for some < ^ < 7r/2. 

We can express G in the basis \a) ,\w): 

G = -U{I -2\0) {0\)U-\I -2 \w){w\) 

= -{I -2U\0){0\U^){I -2\w){w\) 

= -{I-2\a){a\){I-2\w){w\) 

= -/ + 2|o-)(o-| +2|t(;)(t(;| -43;|(T)(t(;| , 
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whence 



G\a) = {I- Ax^)\a) + 2x\w) , 
G\w) = -2x\a) + \w). 

In view of (0), we solve the equation 

smirit) 

— \L1 = 2x = 2 cos 
sm 

for t to get the solution 

TT — 29 11 — 26 TT — 2 arccos x 

as in (|6|). It is then easy to check that e"*''^*" \(t) = G\ct) and that e~*''^*° \ w) = G\w). 

Let S be the subspace spanned by \a) and \w) and let S"-*" be its orthogonal complement. Let 
P be the orthogonal projection onto 5"*". We have just shown that e"*^*" = G restricted to S. For 
I a) G S""*", clearly G|a) = — while e~^^^° leaves pointwise fixed. Thus, 

^-iHto = G + 2P, 

and since GP = PG = -P = -P^, we have e~^*^*o = as expected. 



Remark. By adding f-P to H, we get a slightly more complicated Hamiltonian H such that 



-iHt, 



to 

° = G on the whole Hilbert space. 



Finally, we show how close to is to 1, assuming x << 1. Expanding to as a power series in x, 
we get 



to = l + -x2 + 0(x^), 



and thus 



-iH{2x^/3+0{x'^)) 



--iHx^ + 0{x^)H 
3 



-xVl - 2;2 + O(x^) 
-x3 + 0(x5). 



The second equation comes from expanding the exponential as a power series. The third equation 
holds because \H\ = |^sin20| = xVT^-a?. When x = N~'^^'^, we see that e~*^ comes within 
0(iV-3/2) of G + 2P in norm, and thus e'^'^ comes within 0(iV-3/2) q£ q2 



Corollary 2 If t = f A^^/^ t/ien e"*^* |o-) = ju-) + 0(l/iV). 
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4.3 Discussion 



Farhi and Gutmann show that the Hamiltonian H' finds the target in optimal time in the fohowing 
sense: no other Hamiltonian of the same form — that is, H^ + E \w){w\ where Hd has no special de- 
pendence on w — can find \ w) any faster, even if is allowed to depend on time. Om' Hamiltonian 
is clearly not of this form, so their lower bounds aren't directly applicable here. Indeed, it is only 
by the lower bounds shown for "digital" quantum search that we know that our Hamiltonian is 
optimal to simulate a small digital quantum circuit. It is an interesting question whether one can 
deduce the same lower bound by more direct means. 



5 Further Research and Open Problems 

We have seen how Grover's algorithm can be described much more simply using a Hamiltonian than 
directly with unitary operators. We don't present details them here, but variants of Grover's original 
algorithm also admit simple Hamiltonian descriptions. There may be new, yet unknown quantum 
algorithms which are more easily described with Hamiltonians than with unitary operators, and 
which may indeed be first discovered by their Hamiltonians. 

There are two principal challenges to fashioning a new quantum algorithm via a Hamiltonian: 

1. finding an appropriate, and hopefully intuitive, Hamiltonian for the problem at hand, and 

2. deciding how (or if) the time evolution governed by such a Hamiltonian approximates a true 
(digital) quantum algorithm given by a small quantum circuit. 

In the case of Grover's algorithm considered here, we were fortunate to achieve both goals. Grover's 
algorithm came first, however, so we knew what to shoot for. Even so, the intuition provided 
in Section ^ may be useful for constructing new algorithms, or at least viewing other existing 
algorithms from a different angle. 

Our results thus point to an important general question: when, given a Hamiltonian on a 
system of qubits, can the corresponding time evolution be simulated (even approximately) by a 
small quantum circuit? Is there an easy criterion, based on the structure of the Hamiltonian itself? 
Such a criterion would provide a new way to prototype new quantum algorithms via Hamiltonians. 

Can Farhi &; Gutmann's original H' be simulated efficiently by a quantum circuit? 

Recently, Farhi, et al. show how to solve certain instances of SAT with slowly time-dependent 
Hamiltonians (adiabatic evolution). Their results provide good physical intuition. Is there a corre- 
sponding algorithmic intuition? 

Can one find an intuitive Hamiltonian for a quantum factoring algorithm? 
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